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Abstract. Hyperbolic polynomials are real polynomials whose real hypersurfaces are nested 
ovaloids, the inner most of which is convex. These polynomials appear in many areas of 
mathematics, including optimization, combinatorics and differential equations. Here we in- 
vestigate the special connection between a hyperbolic polynomial and the set of polynomials 
that interlace it. This set of interlacers is a convex cone, which we write as a linear slice 
of the cone of nonnegative polynomials. In particular, this allows us to realize any hyper- 
bolicity cone as a slice of the cone of nonnegative polynomials. Using a sums of squares 
relaxation, we then approximate a hyperbolicity cone by the projection of a spectrahedron. 
A multiaffine example coming from the Vamos matroid shows that this relaxation is not 
always exact. Using this theory, we characterize the real stable multiaffine polynomials that 
have a definite determinantal representation and construct one when it exists. 

I. Introduction 

A homogeneous polynomial / G M[x] of degree d in variables x = (x\,. . . ,x n ) is called 
hyperbolic with respect to a point e G R n if /(e) 7^ and for every a G M. n , all roots of the 
univariate polynomial f(te + a) G are real. Its hyperbolicity cone, denoted C(f,e) is 
the connected component of e in M. n \V^,(f) and can also be defined as 

C(f, e) = {a G M n : f(te - a) ^ when t < 0}. 

As shown in Garding [6], C(f, e) is an open convex cone and / is hyperbolic with respect 
to any point contained in it. Hyperbolicity is reflected in the topology of the real projective 
variety Vr(/) in P n_1 (IR). If Vr(/) is smooth, then / is hyperbolic if and only if Vr(/) 
consists of |_f J nested ovaloids, and a pseudo-hyperplane if d is odd (see [HI Thm. 5.2]). 




FIGURE 1. A quartic hyperbolic hypersurface and two of its affine slices. 

A hyperbolic program, introduced and developed by Giiler [7J, Renegar [H] and others, 
is the problem of maximizing a linear function over an affine section of the convex cone 
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C(f,e). This provides a very general context in which interior point methods are effective. 
For example, taking f = Yli x i an d e = (1, . . . , 1), we see that C(f, e) is the positive orthant 
(]R + ) n and the corresponding hyperbolic program is a linear program. If instead we take / as 



the determinant of a symmetric matrix of variables X 
then C(f, e) is the cone of positive definite matrices. 



and e is the identity matrix, 



/ 


e 


C(f,e) 


hyperbolic program 




(1,...,1) 


(R+)» 


linear program 


det(X) 


I 


positive definite matrices 


semidefinite program 



It is a fundamental open question whether or not every hyperbolic program can be rewrit- 
ten as a semidefinite program. Helton and Vinnikov j8] showed that if / G M[xi, £2, £3] 
is hyperbolic with respect to a point e, then / has a definite determinantal representation 
f = det(^XjMj) where M\, M%, M3 are real symmetric matrices and the matrix e»Mj 
is positive definite. Thus every three dimensional hyperbolicity cone is a slice of the cone 
of positive semidefinite matrices. For a survey of these results and future perspectives, see 
also [18J. On the other hand, Branden [2J has given an example of a hyperbolic polynomial 



/ (see Example 5.11) such that no power of / has a definite determinantal representation. 
There is a close connection between definite determinantal representations of a hyperbolic 
polynomial / and polynomials of degree one-less that interlace it, which has also been used 
in |12j to study Hermitian determinantal representations of hyperbolic curves. 



Definition 1.1. Let f,g G R[t] be univariate polynomials with only real zeros and with 
deg(g) = deg(/) — 1. Let aii ^ • • ■ ^ be the roots of /, and let /?i ^ • • • ^ fid-i be the 
roots of g. We say that g interlaces / if a, ^ 0i ^ cti+i for all % — 1, . . . , d — 1. If all these 
inequalities are strict, we say that g strictly interlaces /. 

If / G M.[x] is hyperbolic with respect to e and g is homogeneous of degree deg(/) — 1, we 
say that g interlaces / with respect to e if g(te + a) interlaces f(te + a) for every a G M. n . 
This implies that g is also hyperbolic with respect to e. We say that g strictly interlaces 
/ if gite + a) strictly interlaces f(te + a) for a in a nonempty Zariski-open subset of IR n . 

The most natural example of an interlacing polynomial is the derivative. If f(t) is a real 
polynomial with only real roots, then its derivative f'(t) has only real roots, which interlace 
the roots of /. Extending this to multivariate polynomials, we see that the roots of ^_f(te+a) 
interlace those of f(te + a) for all a G R n . Thus 

n 

Def = £' 



df_ 

'dxi 



interlaces / with respect to e. If / is square-free, then D e f strictly interlaces /. This was 
already noted by Garding [6 J and has been used extensively, for example in [1| and [14] : for 
general information on interlacing polynomials, see also [5J and [13, Ch. 6]. 

Remark 1.2. If / is square-free and d = deg(/), then / is hyperbolic with respect to e if 
and only if f(te + a) has d distinct real roots for a in a Zariski-open subset of M. n . In this 
case, if g interlaces / and has no common factors with /, then g strictly interlaces /. 
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FIGURE 2. Two affine slices of a cubic interlacing a quartic. 

In this paper, we examine the set of polynomials in R[x]d-i interlacing a fixed hyperbolic 
polynomial. The main result is a description of a hyperbolicity cone C(f, e) as a linear slice 
of the cone of nonnegative polynomials. Using the cone of sums of squares instead gives an 
inner approximation of C(f, e) by a projection of a spectrahedron. This is closely related to 
recent results due to Netzer and Sanyal [10] and Parrilo and Saunderson [11]. We discuss 
both this theorem and the resulting approximation in Section [3} In Section [4] we see that the 
relaxation we obtain is exact if some power of / has a definite determinantal representation. 
A multiaffme example for which our relaxation is not exact is discussed in Section [5] Here we 
also provide a criterion to test whether or not a hyperbolic multiaffme polynomial has a defi- 
nite determinantal representation. The full cone of interlacers has a nice structure, which we 
discuss in Section |6j First we need to build up some basic facts about interlacing polynomials. 

Acknowledgements. We would like to thank Alexander Barvinok, Petter Branden, 
Tim Netzer, Rainer Sinn, and Victor Vinnikov for helpful discussions on the subject of this 
paper. Daniel Plaumann was partially supported by a Feodor Lynen return fellowship of 
the Alexander von Humboldt-Foundation. Cynthia Vinzant was partially supported by the 
National Science Foundation RTG grant DMS-0943832 and award DMS-1204447. 

2. Interlacers 

Let / be a homogeneous polynomial of degree d that is hyperbolic with respect to the 
point e G W 1 . We will always assume that /(e) > 0. Define Int(/, e) to be the set of real 
polynomials of degree d — 1 that interlace / with respect to e and are positive at e: 

Int(/, e) = {g G M[rr]d_i : g interlaces / with respect to e and g(e) > 0}. 

As noted above, the hyperbolicity cone C(f, e) depends only on / and the connected compo- 
nent of W l \ Vr(/) containing e. In other words, we have C(f, e) = C(f, a) for all a G C(f,e). 
We will see shortly that Int(/, e) does not depend on e either, but only on C(f, e). 

Theorem 2.1. Let f G M.[x]d be square-free and hyperbolic with respect to e G lR n , where 
/(e) > 0. For h G IRfx^-i, the following are equivalent: 

(1) helnt(f,e); 

(2) h G Int(/,a) for all a G C(f,e); 

(3) D e f ■ h is nonnegative on Vr(/); 

(4) D e f ■ h — f ■ D e h is nonnegative on R n . 
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The proof of this theorem and an important corollary are at the end of this section. First, 
we need to build up some theory about the forms in Int(/, e). 

Lemma 2.2. Suppose f\, f 2 , and h are real homogeneous polynomials. 

a) The product f\ ■ f 2 is hyperbolic with respect to e if and only if both fi and f 2 are hyperbolic 
with respect to e. In this case, C(f\ ■ / 2 , e) = C(fi, e) fl C(f 2 , e). 

b) If fi and / 2 are hyperbolic with respect to e, then fi ■ h interlaces fx ■ f 2 if and only if h 
interlaces f 2 . 

c) If h interlaces (fi) k f 2 for fc£N, then (/i) fe_1 divides h. 

Proof. These statements are checked directly after reducing to the one-dimensional case. □ 

Lemma 2.3. For any g and h in Int(/, e), the product g ■ h is nonnegative on Vr(/). 

Proof. To prove this statement, it suffices to restrict to any line x = te + a where a G M. n . 
Suppose that f(te + a) G M.[t] has roots oc\ < . . . < ota and g(te + a) and h(te + a) have roots 
Al < • • • < (3d- 1 and 71 < . . . < 7d-i, respectively. By the assumption that both g and h 
interlace /, we know that 7, G [ai, c^+i] for all 1 < i < g? — 1. Thus, if aj and atj are not 
also roots of g(te + a) or h(te + a), the polynomial gite + a)h(te + a) has an even number of 
roots in the interval [aii, ctj}. Then the sign of giptie + a)h{otie + a) is the same for all % for 
which it is not zero. Because g(e)h(e) > 0, we see that that sign must be nonnegative. □ 




FIGURE 3. Affine slices of two cubics interlacing a quartic. 



Lemma 2.4. Suppose that f is square-free and that g G Int(/, e) strictly interlaces f. Then 
a polynomial h G Rfarj^-i belongs to Int(/, e) if and only if g ■ h is nonnegative on Vr(/). 



Proof. One direction follows from Lemma 2.3 For the other, let h G IR[x]d-i for which g ■ h 
is nonnegative on Vr(/). First, let us consider the case where / and h have no common 
factor. Then, for generic a G M n , the roots of f(te + a) are distinct from each other and from 
the roots of g(te + a) and h(te + a). The product g(te + a)h(te + a) is then positive on all 
of the roots of fite + a). Since g(te + a) changes sign on consecutive roots of fite + a), we 
see that h(te + a) must have a root between each pair of consecutive roots of f(te + a), and 
thus h interlaces / with respect to e. 

Now suppose / = fi ■ f 2 and h — fi ■ h±. We will show that hi interlaces f 2 , and thus h 
interlaces /. Again, we can choose generic a for which the roots of f(te + a) and g(te + a) 
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V(f) 



FIGURE 4. Non-crossing arcs of Lemma 2.6 



are all distinct. Consider two consecutive roots a < of the polynomial f'2{te + a). Let 
A; be the number of roots of fi(te + a) in the interval (a,/3). Because g strictly interlaces 
f = fi ■ /2, its restriction g>(te + a) must have k + 1 roots in the interval (a, Thus the 
polynomial g(te + a)fi(te + a) has an odd number of roots in this interval and must therefore 
have different signs in a and (3. Since g ■ fi ■ h\ > on V(/), the polynomial 7ii(ie + a) must 
have a root in this interval. Thus h\ interlaces f 2 and h interlaces /. □ 

Example 2.5. In the above lemma, it is indeed necessary that / and g be without common 
factors. For example, consider / = (x 2 + y 2 — z 2 )(x — 2z) and g = (x 2 + y 2 — z 2 ). Both / 
and g are hyperbolic with respect to [0 : : 1] and g interlaces / with respect to this point. 
However if h = y(x — 2z), then g ■ h vanishes identically on Vr(/) but h does not interlace /. 

Lemma 2.6. For a £ C(f,e), the polynomial D e f ■ D a f is nonnegative on Vr(/). 

Proof. For any b £ C(f,e) and x £ Vr(/), let ai(b,x) < ... < aa(b,x) denote the roots 
of f(tb + x). Because C(f,e) is convex, the line segment joining e and a belongs to this 
cone. As we vary b from e to a along this line segment, the roots {a:j(&, x)b + x}j g ^], form d 
non-crossing arcs in the plane x + span{e, a}, as shown in Figure |4j Since fix) = 0, one of 
these arcs is just the point x. That is, there is some k for which x) = for all b in the 
convex hull of e and a. 

Now /(e) > implies f[b) > for all b £ C(f,e). Thus §j.f{tb + sc) is positive for 
t > ad(b,x). Furthermore, the sign of this derivative on the ith root, ati(b,x) depends only 
on i. Specifically, for alH = 1, . . . , d, 




□ 



Now we are ready to prove Theorem 2.1 
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Proof of Theorem 2.1 , 
(4 =^ 3) Clear. 

(1 3) If / is square free, then D e f strictly interlaces /. This equivalence then follows 
from Lemma [2.41 

(1, 3 4) Here we need a useful fact about the Wronskian of univariate polynomials. The 
Wronskian of two univariate polynomials p(t), q(t) is the polynomial 



W(p, q) = V ■ q' — p' • q = q 



Hi 

It is a classical fact that if the roots of p and q are all distinct and interlace, then W(p, q) is 
is a nonnegative or nonpositive polynomial [T9| §2.3]. Thus if h G Int(/, e) is coprime to /, 
then for generic x, the roots of f(te + x) and h(te + x) interlace and are all distinct. Thus 
their Wronskian h(te + x)f'(te + x) — h'(te + x)f(te + x) is either nonnegative or nonpositive 
for all t. By (3), the product h(te + x)f'(te + x) is nonnegative on the zeroes of /, so we 
see that the Wronskian is nonnegative. Setting t = gives us that h ■ D e f — D e h ■ f is 
nonnegative for all x G M. n , as desired. If / and h share a factor, say / = fi ■ f'2, h = fx ■ hi, 
we can use the identity W(fi ■ f'2, fi • hi) = ffW(f2, hi) to reduce to the coprime case. 

(2 <^ 1) Because / is square free, both D e f and D a f share no factors with /. Thus D e f 
strictly interlaces / with respect to e and D a f strictly interlaces / with respect to a. 



Suppose h interlaces / with respect to a and h(a) > 0. By Lemma 2.4, h ■ D a f is 



nonnegative on Vr(/). Using Lemma 2.6, we see that D e f ■ D a f is also nonnegative on 
Vr(/). Taking the product, it follows that (D a f) 2 ■ D e f ■ h is nonnegative on Vr(/). Because 
D a f and / have no common factors, we can conclude that D e f ■ h is nonnegative on Vr(/). 
Using Lemma 2.4 again, we have h G Int(/, e). Switching the roles of a and e in this argument 
gives the reverse implication. □ 

Corollary 2.7. The set Int(/, e) is a closed convex cone. If f is square-free, this cone is 
linearly isomorphic to a section of the cone of nonnegative polynomials of degree 2deg(/) — 2: 

(2.1) Int(/,e) = {h G R[a;]deg(/)-i : DJ-h-f-D e h >0 onR n }. 

If f — fi ' f2 where V(f) = V{f2) and f 2 is square-free, then 

Int(/,e) = A -IntC/^e) 

and is isomorphic to a section of the cone of nonnegative polynomials of degree 2 deg(/2) — 2. 



Proof. For square-free /, the description (2.1) follows directly from Theorem 2.1 The map 

h ^ DJ -h- f ■ D e h 

is a linear map from ]R[x]d cg (/)-i to R[x](2deg(/)-2)- We see that Int(/, e) is the preimage of 
the cone of nonnegative polynomials in M[a;](2deg(/)-2) under this map. We can also check 
that this map is injective. Because / is square free, D e f and / are coprime. Hence if / were 
to divide D e f ■ h, then / would have to divide h, which it cannot. Thus D e f • h — f ■ D e h 
cannot be identically zero. 

If / is not square-free, then / factors as fi ■ f 2 as above. By Lemma 2.2 c), any polynomial 
that interlaces / must be divisible by f\. By part (b), the remainder must interlace f 2 . Thus 
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Int(/, e) C fi ■ Int(/2, e). Similarly, if h interlaces f'2, then f\ ■ h interlaces / = f\ ■ fi- Thus 
Int(/, e) is the image of the convex cone Int(/2, e) under a linear map, namely multiplication 
by fx. This shows that it is linearly isomorphic to a section of the cone of nonnegative 
polynomials of degree 2deg(/ 2 ) — 2. □ 

3. HYPERBOLICITY CONES AND NONNEGATIVE POLYNOMIALS 

An interesting consequence of the results in the preceding section is that we can recover 
the hyperbolicity cone C(f, e) as a linear section of Int(/, e), and thus as a linear section of 
the cone of nonnegative polynomials. We show this by considering which partial derivatives 



D a (f) interlace /. Using Theorem 2.1, we often come across the polynomials 

A e ,J = DJ.DJ-f-D e DJ. 

Notice that A e a / is homogeneous of degree 2d — 2, symmetric in e and a, and linear in each. 

Theorem 3.1. Let f G Rfscjd be square-free and hyperbolic with respect to the point e G M. n . 
The intersection of Int(/, e) with the plane spanned by the partial derivatives of f is the 
image of C(f, e) under the linear map a 1— > D a f . That is, 



(3.1) C(f,e) = {aeR n : DJ G Int(/,e)}. 



Furthermore, C(f, e) is a section of the cone of nonnegative polynomials of degree 2d — 2: 



(3.2) C(f,e) = {aeR n : A e ,J > on R n }. 



Proof. Let C be the set on the right hand side of (3.1 ). From Theorem 2.1 we see that D a f 
interlaces / with respect to e for all a G C(f,e). This shows C(f,e) C C and hence the 
inclusion C(f,e) C C, since C is closed. If this inclusion were strict, there woul d exist a 
point a G C\C(f,e) with f(a) ^ 0, since C is also a convex cone by Corollary 



2.7 



Thus to 

show the reverse inclusion, it therefore suffices to show that for any point a outside of C(f, e) 
with f{a) 7^ 0, the polynomial D a f does not belong to Int(/, e). If a belongs to —C(f,e), 
then —D a f belongs to Int(/, e). In particular, —D a f(e) > and D a f does not belong to 
Int(/, e). Thus we may assume a ^ C(f,e) U —C(f,e). Since / is hyperbolic with respect 
to e, all of the roots ati < ... < of f(te + a) are real. The reciprocals of these roots, 
1/aii, . . . 1/otd, are roots of the polynomial f(e + ta). 

Because a is not in C(f, e) U —C(f, e), there is some 1 < i < n for which cnj < < a, l+ i. 
Since /(e) 7^ and f(a) 7^ 0, we can take reciprocals to find the roots of f(e + ta): 

11 111 1 
— < <...< — <0< — < < ... < . 



By Rolle's Theorem, the roots of J^/(e + ta) interlace those of /(e + ta). Note that the 
polynomial §~ t f(e + ta) is precisely D a f(e + ta), so the roots of D a f(e + ta) interlace those of 
f(e + ta). In particular, there is some root /3 of D a f(e + ta) in the open interval (1/cci, l/ad), 
and thus 1//3 G" [a\, a d \ is a zero of D a f(te + a). Therefore D a f(te + a) has only d — 2 roots 
in the interval [ct\, a.j\ and thus cannot interlace / with respect to e. 



Combining this with Theorem 2.1 shows the equality in (3.2). □ 
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FIGURE 5. For a outside of the hyperbolicity cone, D a f does not interlace /. 



Corollary 3.2. Relaxing nonnegativity to sums- of- squares in (3.2) gives an inner approxi- 
mation to the hyperbolicity cone of f : 

(3.3) {oSR" : A eta f is a sum of squares } C C(f,e). □ 



If the relaxation (3.3) is exact, then the hyperbolicity cone is a projection of a spectra- 



hedron, namely of a section of the cone of positive semi definite matrices in , where 

N = ( n+ _7 2 ) = dimR[x](2_i. A polynomial F is a sum of squares if and only if there exists 
a positive semi definite matrix G such that F = v T Gv, where v is the vector of monomials 
of degree at most deg(F)/2. We call such a matrix G a Gram matrix of F. The linear 
equations giving the Gram matrices of A e ^ a f give the desired section of IR^^. 



If the relaxation (3.3) is not exact, one can allow for denominators in the sums of squares 



and successively improve the relaxation. More precisely, for any integer N ^ consider 
jflGR" : xf) N ■ A e ^f is a sum of squares } C C(f,e). 



(3.4) 



C 



N 



As above, Cjv is seen to be a projection of a spectrahedron. Furthermore, by a result of 
Reznick in [J5], for any positive definite form F £ M.[x] there exists some positive integer N 
such that (X^ILi x i) ■ F is a sum of squares. Thus if Vr(/) is smooth, then {A e a f \ a £ M. n } 
contains a strictly positive polynomial, for example A e>e f. It follows that the hyperbolicity 
cone C(f, e) is the closure of the union of all the cones Cn- 

Remark 3.3. In a recent paper |10j . Netzer and Sanyal showed that the hyperbolicity cone 
of a hyperbolic polynomial without real singularities is the projection of a spectrahedron. 
Their proof uses general results on projected spectrahedra due to Helton and Nie and is not 



fully constructive. In particular, it does not imply anything about equality in (3.3) or (3.4). 



Explicit representations of hyperbolicity cones as projected spectrahedra have recently 
been obtained by Parrilo and Saunderson in [11] for elementary symmetric polynomials and 
for directional derivatives of polynomials possessing a definite determinantal representation. 

Remark 3.4. We also have the relaxation 

(3.5) {a £ IR n : D e f ■ D a f is a sum of squares modulo (/)} C C(f,e). 



It is unclear whether or not this relaxation is always equal to (3.3). Its exactness would 



also show C(f, e) to be the projection of a spectrahedron. We will see below that if / has a 



definite determinantal representation, then we get equality in (3.3) and (3.5) 
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Example 3.5. Consider the quadratic form f(x) 



x^, which is hyperbolic 



with respect to the point e = (1, 0, . . . , 0). The hyperbolicity cone C(f, e) is known as the 
Lorentz cone. In this example, the relaxation ( |3.3[ ) is exact. To see this, note that 

(2xi)(2a!a;i — s ^^2ajXj) — {x\ — x 2 ){2ai) 
2(aix\ — 2 ctjXiXj + aix 2 ). 



Since every nonnegative quadratic form is a sum of squares, there is equality in (3.3). In 
fact, taking the Gram matrix of |A eja /, we recover the Lorentz cone as 

/ a\ —ci2 . . . — a n \ 
— d2 ai 



C(f, e) 



< 



a e 



-a. 



y o 



> 



... ai ) 

Note also that this Gram matrix gives a definite determinantal representation of a\~ 2 f{a). 
Example 3.6. Consider the hyperbolic cubic polynomial 

f = (x- y)(x + y)(x + 2y) - xz 2 } 

with e = [1 : : 0]. Here the polynomial A eja / has degree four in x,y,z. In this case, 
the relaxation (3.3) is exact, as shown in Corollary 4.5 below. (One can also see exactness 
from the fact that every nonnegative ternary quartic is a sum of squares). Using the basis 
(x 2 ,y 2 , z 2 ,xy,xz,yz) of M.[x,y, z]2, we can then write the cone C(f,e) as the set of (a,b,c) 
in IR 3 for which there exists (gi, . . . , £ ^ 6 to make the real symmetric matrix 



/3a + 2b 


9i 


92 


4a - 26 


-2c 


93 \ 


9i 


9a + 26 


9i 


4a - 86 


95 


-2c 


92 


94 


a 


96 








4a - 2b 


4a -86 


96 


8a - 206 - 2#i 


—2c — $3 




-2c 


95 





-2c - g 3 


26 - 2g 2 


-2a - g 6 


\ 93 


-2c 





-95 


-2a - g 6 


2a + 66 - 2gJ 



positive semidefinite. 



The sums of squares relaxation (3.3) is not always exact. A counterexample comes from 
a multilinear hyperbolic polynomial and will be discussed in Example |5.11| 



4. Definite Symmetric Determinants 

We consider det(X) as a polynomial in M[Xjj : i < j G [d]], where X = (Xij) is a 
symmetric matrix of variables. Since all eigenvalues of a real symmetric matrix are real, 
this polynomial is hyperbolic with respect to the identity matrix. The hyperbolicity cone 
C(det(X),J) is the cone of positive definite matrices. Hence, for any positive semidefinite 
matrix E ^ 0, the polynomial 

(4.1) D E {det{X)) = tr(£-X adj ) 
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interlaces det(X). This holds true when we restrict to linear subspaces. For real symmetric 
d x d matrices Mi, . . . , M n and variables x = (xi, 



,x n ), denote 



M(x) 



If M(e) is positive definite for some e G IR n , then the polynomial det(M(x)) is hyperbolic 
with respect to the point e. 

Proposition 4.1. If M is a real symmetric matrix of linear forms such that M(e) y 
for some e G M. n , then for any positive semidefinite matrix E , the polynomial tr(E ■ M ad j) 
interlaces det(M) with respect to e. 

Proof. By the discussion above, the polynomial D^(det(X)) — tr(E- X ad ->) interlaces det(X) 
with respect to E. (In fact these are all of the interlacers of det(X). See Example 6.2 below.) 
By Theorem 2.1, tr(E • X^) interlaces det(X) with respect to any positive definite matrix, 
in particular M(e). Restricting to the linear space {M(x) :xeR n } shows that tr(E ■ M adj ) 
interlaces det(M) with respect to e. □ 

Theorem 4.2. If f G IR[x]d has a definite symmetric determinantal representation f = 
det(M) with M(e) y and M(a) y 0, then A e a / is a sum of squares. In particular, there 
is equality in (3.3). 

Proof. Because M(e) and M(a) are positive semidefinite, we can write them as sums of 
rank-one matrices: M(e) = £\ AjAf and M(a) = Ylj^j^Ji where Xi,fij G M d . Then 
DJ = (M(e), M ad J) = (£. KXJ, M^) = XjM^\, so 

DJ = XjM^Xi and, similarly, DJ = ^ /'j 



Furthermore, by Proposition 4.6 below, the second derivative D a D^f is 



D e DJ =D e ( J2^J MaA ^J ) = ^2 U H where U H 



M A, y,j 
Xj 
nj 



Now, again using Proposition |4.6| we see that A ejQ / equals 
(4.2) 



J2( (XjM^X^jM^^) - det(M) ■ 



U; 



which is the desired sum of squares. 



□ 



In fact, something stronger is true. We can also consider the case where some power of / 
has a definite determinantal representation. This is particularly interesting because taking 
powers of a hyperbolic polynomial does not change the hyperbolicity cone. 

Corollary 4.3. If f G IR[x]<j and a power f r has a definite symmetric determinantal repre- 
sentation f r = det(M) with M(e), M(a) y 0, then A e>a (/) is a sum of squares. In particular, 
there is equality in (3.3). 



HYPERBOLIC POLYNOMIALS, INTERLACERS, AND SUMS OF SQUARES 11 

Proof. Let f r have a definite determinantal representation. We have A eia (/' r ) = r/ 2 ^ -1 ^ A eo /. 
Theorem 



4.2 



states that A e;a (/ r ) is a sum of squares, 

9! + ■ ■ • + 9 2 s = rf^A e J 

for some gi G IR[x]. Let p be an irreducible factor of J 2 ^" 1 ). Then p is hyperbolic with 
respect to e and the right hand side vanishes on V<c{p). Therefore, each gi vanishes on Vr(p) 
and thus on Vc(p), since V^(p) is Zariski dense in Vc{p)- Thus we can divide the g^ by p. 
By iterating this argument, we get the claim. □ 

Remark 4.4. This result is closely related to (but does not seem to follow from) [9j 
Thm. 1.6], which says that the parametrized Hermite matrix of / is a sum of matrix squares 
whenever a power of / possesses a definite determinantal representation. 



Corollary 4.5. If f £ M.[xi, X2, £3], then there is equality in (3.3). 

Proof. By the Helton- Vinnikov Theorem [8], every hyperbolic polynomial in three variables 



has a definite determinantal representation. The claim then follows from Theorem 4.2 □ 



The following determinantal identities were needed in the proof of Theorem 4.2 above. 
Proposition 4.6. Let X be a d x d matrix of variables Xij and let \ ■ \ denote det(-). Then 



for any vectors a 


P,7,8 


e 


C d 


we 


have 




(4.3) 


X f3 




X 


5 




X 5 




X /3 


a T 




T 

7 







a T 




7 T 



IXI 



X (3 5 
T 




a 



in C[Xij : 1 < i, j < d\. Furthermore, 



D/3a T \X\ 



X 

T 

a 



and D Sl TDp a T\X 



7 



X 13 5 
a T 
7 T 



Proof. We will prove the first identity using Schur complements. See, for example, [3, §1]. 



If A is a m x m submatrix of the n x n matrix 



B D 



\A\ ■ \D — BA 1 C\. If D is the zero matrix, this simplifies to 



A C 
B 



\A\ 



-1 

L4J 



• BA^C 



\A\ 



To obtain the desired identity, we set A = X, B = 

^X adj (/3 5) 



then its determinant equals 

■ \BA adi C\. 
and C ={13 5): 



X (3 S 
a T 
T 



7 



\X\ 



X 



a 

1 T ) 



1 

\X\ 



a T X^(3 a T X^5 



Multiplying both sides by det(X) finishes the proof of the determinantal identity. 

For the claim about derivatives of the determinant, by additivity, we only need to look at 
the case when a, (3, 7, 5 are unit vectors, e,, e^, e^, ei, respectively. Then Dp a T\X\ = D e . e r\X\ 
is the derivative of \X\ with respect to the entry X^. This is the signed minor of X obtained 
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by removing the jth row and ith column, which is precisely the determinant 
the derivative of this determinant with respect to Xik the same way gives 



X ej 
eT 



<9 2 |X| 
dXjidXik 



D^ T D e 



\x\ 



X 







X ej ei 
ef 
el 



. Taking 



□ 



We conclude this section with a general result inspired by Dixon's construction of determi- 
nantal representations of plane curves, which will be applied in the next section. If / G R[x]d 
has a definite determinantal representation, / = det(M) with M(e) >~ 0, then M ad J is a d x d 
matrix with entries of degree d — 1. This matrix has ran k at most one on V(f), as seen by 
the identity M ■ M adj = det(M) • /. By Proposition 4.1, the top left entry interlaces 
/ with respect to e. In fact, these properties of M ad -> are enough to reconstruct a definite 
determinantal representation M. 



j) be a symmetric d x d matrix of real forms of degree d — 1. 



1 . If A has rank 
[l/f d - 2 )A^ has 



Theorem 4.7. Let A = 

Suppose that f G is irreducible and hyperbolic with respect to e G 

one modulo (/), then f d ~ 2 divides the entries of A ad i and the matrix M ■ 
linear entries. Furthermore there exists 7 G ffi such that 

det(M) = 7/. 

If an interlaces f with respect to e and there is no vector v G IR n \{0} in the kernel of A, 
then 7 7^ and M(e) is definite. 

Proof. By assumption, / divides all the 2x2 minors of A. Therefore, f d ~ 2 divides all of 
the (d — 1) x (d — 1) minors of A and thus all of the entries of the adjugate matrix A^, 
see [12| Lemma 4.7]. We can then consider the matrix M = (l/f d ~ 2 ) ■ A ad K By similar 
arguments, f d ~ x divides det(A). Because these both have degree d(d — 1), we conclude that 
det(.A) = A/ d_1 for some A G R. Putting all of this together, we find that 

det W = 7^y det (^ adJ ) = 7^ det (^ = ^ dl f ' 

so we can take 7 = A^" 1 . Now, suppose that an interlaces / and that 7 = A = 0. Then 
det(A) is identically zero. In particular, the determinant of A(e) is zero, there is some 
nonzero vector v G IR n in its kernel, and v T A(e)v is also zero. 

We will show that the polynomial v T Av is not identically zero and that it interlaces / 
with respect to e. This will contradict the conclusion that v T A(e)v = 0. Because A has rank 
one on V(/), for each i — 1, . . . , d we have that 

(4.4) (ejA ei )(v T Av)-(ejAv) 2 = modulo (/). 

If v T Av is identically zero, then ejAv vanishes on V(/). Since ejAv only has degree d— 1, it 
must vanish identically as well. As this holds for each i, this implies that Av is zero, which 
contradicts our assumption. Thus v T Av cannot be identically zero. 

Furthermore, (4.4) shows that an • (v T Av) is nonnegative on Vr(/). Then Lemma 2.4 



shows that v T Av interlaces / with respect to e. In particular, v T Av cannot vanish at the 
point e. Thus the determinant of A and hence M cannot be identically zero. 
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Thus M is a determinantal representation of /. To show that M(e) is definite, it suffices 
to show that A(e) is definite. For any vector v G K™, we see from ( 4.4[ ) with % — 1 that 
a\\V T Av is nonnegative on Vr(/). Thus an(e) -t^Au belongs to Int(/, e) by Lemma 2.4 and 
in particular an(e) ■ f T A(e)i> is positive for all v G M d . Hence the matrix A(e) is definite. □ 



5. MULTIAFFINE POLYNOMIALS 

An interesting special case of a hyperbolic polynomial is a multiaffine polynomial whose 
hyperbolicity cone contains the positive orthant. These polynomials are deeply connected 
to the theory of matroids [2| 14] [20]. 

Definition 5.1. A polynomial / G lR[x] is called affine in x^ if the degree of / in Xj is at 
most one. If / is affine in each variable x\, . . . , x n , then / is called multiaffine. 

Much of the literature on these polynomials deals with complex polynomials, rather than 
real polynomials, and the property of stability in place of hyperbolicity. 

Definition 5.2. A polynomial / G C[x] is called stable if f(fi) is non-zero whenever the 
imaginary part of each coordinate /ij is positive for all 1 < % < n. 

A real homogeneous polynomial / G M[x] is stable if and only if / is hyperbolic with 
respect to every point in the positive orthant. After a linear change of variables, every 
hyperbolic polynomial is stable. In 2004, Choe, Oxley, Sokal, and Wagner [3] showed that if 
/ G is stable, homogeneous, and multiaffine, then its support (the collection of / C [n] 

for which the monomial Yliei Xi a PP ears m /) is the set of bases of a matroid. They further 
show that any representable matroid is the support of some stable multiaffine polynomial. 
In 2010, Branden [2] used this deep connection to disprove the generalized Lax conjecture 



by showing that the bases generating polynomial of the Vamos matroid (see Example 5.11) 



is hyperbolic but none of its powers has a determinantal representation. This example will 



also provide a counterexample to equality in our relaxation (3.3). 

The Wronksian polynomials A e>a / also played a large part in the study of multiaffine 
stable polynomials. They are particularly useful when the points e and a are unit vectors. 
In this case, we will simplify our notation and write 

df df d 2 f 
Ay(/) := X.,M) = i),./ i)X] - ■ i)Xii)x - 

Using these polynomials, Branden [lj established a necessary and sufficient condition for 
multiaffine polynomials to be stable. 

Theorem 5.3 ([T], Theorem 5.6). For multiaffine f G M.[x], the following are equivalent: 

(1) Ay/ is nonnegative on W 1 for all 1 <i,j <n, 

(2) f is stable. 

Branden also notes that the implication (2)=>(1) holds for polynomials that are not mul- 
tiaffine, giving an alternative proof of a part of Theorem |2.1| above. The other implication 



however, does not hold in general, as the following example shows. 
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Example 5.4. Let h = x\ + x\, q = X\ + X2 and iVeN. Clearly g /t is not hyperbolic with 
respect to any e G M, 2 , but for all 2, j G {1, 2} we have 

^M N h) = q 2N ^h + Nq 2N - 2 h 2 ^q 

= q 2N ~ 2 (q 2 A l3 h + Nh 2 ). 

Now let z G R be the minimal value that q 2 Aijh takes on the unit sphere and let N > \z\. 
Then, since Aij(q N h) is homogeneous, we see that Aij(q N h) is nonnegative on Mr. Because 
Aij(q N h) is a homogeneous polynomial in two variables, it is even a sum of squares. Thus 
A a b{q N h) is a sum of squares for all a, b in the positive orthant. This also shows that the 



converse of Corollary 4.3 is not true, i.e. there is some polynomial / such that A ea f is a 
sum of squares for all e, a in some full dimensional cone, but no power of / has a definite 
determinantal representation. 

In an analogous statement, the polynomials Ajj can also be used to determine whether or 
not a homogeneous multiaffme stable polynomial has a definite determinantal representation. 

Theorem 5.5. Let f G M.[x]d be homogeneous and stable. Suppose f is affine in the variables 
Xi, . . . , Xd and the coefficient of X\ - ■ ■ Xd in f is non-zero. Then the following are equivalent: 

(1) Aijf is a square in M[x] for all 1 < i, j < d; 

(2) §t . dJL i s a S q U are in R[x]/(f) for all 1 < i,j < d; 

(3) f has a definite determinantal representation. 

Lemma 5.6. Let f G M[x] be affine in Xi and Xj for some i,j G [n]. If f = g ■ h with 
g, h G M[x], then A^f is a square if and only if A^g and A^h are squares. 

Proof. Suppose Ay/ is a square. Since / is affine in Xi,Xj, both g and h are affine in Xi,Xj 
and either = or = 0. It follows that either Aijg = or Aijh = 0. Using the identity 
Aijf = g 2 Aijh + h 2 Aijg, we see that either A^g = or A^g = (Aijf)/h 2 . In both cases A^g 
is a square. The same holds true for A^h. For the converse, suppose that A^g and A^h are 
squares. As we saw above, one of them is zero. Thus Aj,/ = h 2 Aijg or A^f = g 2 A^h. □ 



Proof of Thm. \5^\ (1^2) Clear. 

(2 ^> 3) For a start, supp ose t hat / is irreducible. We will construct a matrix A satisfying 
the hypotheses of Theorem 4.7 For every % < j, the polynomial ■ is equivalent to a 

square ajj modulo (/). In the case i = j we can choose an = J^. Then easy to check that 
anOu equals a 2 Xi modulo (/). Further, for every 2 < i < j < d, the polynomials (and^) 2 
and (auaij) 2 are equivalent modulo /. After changing the sign of a^ if necessary, we see 
that a\\Oij equals auaij modulo (/). Because / is irreducible, it follows that the symmetric 
matrix A = (aij)ij has rank one on V(/). 

We now need to show that there is no vector v G IR n \{0} in the kernel of A. For each 
k = 1, . . . , d, consider the point pk = Ylj£{d]\{k} e h "which hes in the real variety of /. For 
j ^ k, we see that df/dxj vanishes at pk, and therefore so must a^j. On the other hand, 
Ofcfc(Pfc) = df /dxk(pk) equals the nonzero coefficient of X\---Xd in /. Now suppose that 
Av = for some v G M. n . The kth row of this is ^ . VjOkj = 0. Plugging in the point pk then 
shows that must be zero, and thus v is the zero vector. Since / is stable, an = df/dxi 



interlaces it, and so by Theorem 4.7 / has a definite determinantal representation. 
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If / is reducible and g is an irreducible factor of /, then, by Lemma 5.6, A^g is a square. 
Since every irreducible factor of / has a definite determinantal representation, so has /. 

(3 =>- 1) Let / = det(M) = det(^XjMj) where Mi, . . . ,M n are real symmetric d x d 
matrices where £V Mj >- 0. Because / is affine in each of the variables Xi,...,Xd, the 
matrices Mi,..., Ma must have rank one. Furthermore, since / is stable, these rank-one 
matrices must be positive semidefinite (see j2], proof of Theorem 2.2). Thus we can write 



Mi = Vivf, with Vi G M d for each 1 < i < d. Then by (4.2) and Proposition 4.6 we have 
A {j f = (vfM^hj) 2 for 1 < i,j < d. □ 

Corollary 5.7. Let f G WL[x] be homogeneous, stable and multiaffine. Then the following 
are equivalent: 

(1) Ay / is a square for all 1 < i,j < n; 

(2) §^ ■ M~ is a square in M[x]/(f) for all 1 < i,j < n; 

(3) f has a definite determinantal representation. 

Proof. This is an immediate consequence of the preceding theorem. □ 

Corollary 5.8. Let 1 < k < n and let f G M[x] be a multiaffine stable polynomial. If f has a 
definite determinantal representation, then J^- and f\ Xk =o also have a definite determinantal 
representation. 

Proof. Let 1 < k,i,j < n, g — and h = f\ Xk =o- Wagner and Wei j20] calculated 

^ijf = x\ ■ A {j g + x k -p + A i:j h, 

where p,g,h G . . . , x n ] do not depend on x^. Since Ay/ is a square, A^g and Ay/i are 

squares as well. Thus g and h have a definite determinantal representation. □ 

Corollary 5.9. Let f = g ■ h, where f,g,h£ M[x] are multiaffine stable polynomials. Then 
f has a definite determinantal representation if and only if both g and h do. 



Proof. This follows directly from Lemma 5.6| and Theorem 5.5 □ 



Example 5.10 (Elementary Symmetric Polynomials). Let e<2 G M[xi, . . . , x n ] be the ele- 
mentary symmetric polynomial of degree d. We have Ay-ei = 1, A^e n = and Ay-e n _i = 
(xi . . .x n /xiXj) 2 for all 1 < % < j < n. It is a classic result, that these are the only cases 
where has a definite determinantal representation |T6j. Indeed, for n > 4 and 2 < d < n — 2 
the coefficients of the monomials (£3X5 • • • Xd+2) 2 , (^4^5 • • • Xd+2) 2 and ^3X4(^5 • • • Xd+2) 2 in 
A\2ed are all 1. Specializing to Xj = 1 for j > 5 then shows that A12/ is not a square. 



Example 5.11 (The Vamos Polynomial). The relaxation (3.3) is not always exact. An 
example of this comes from the multiaffine quartic polynomial in M.[xi, . . . , Xgk given as the 
bases generating polynomial of the Vamos matroid: 

h(xi,...,x 8 ) = ^2 Yl x ^ 
where C = {{1, 2, 3, 4}, {1, 2, 5, 6}, {1, 2, 7, 8}, {3, 4, 5, 6}, {3, 4, 7, 8}}. Wagner and Wei |2U] 



have shown that the polynomial h is stable, using an improved version of Theorem 5.3 



and representing Ai 3 h as a sum of squares. But it turns out that A 78 h is not a sum of 
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squares. Because the cone of sums of squares is closed, it follows that for some a, e in the 
hyperbolicity cone of h, the polynomial D e h ■ D a h — h ■ D e D a h is not a sum of squares. 
In order to show that /S.^h is not a sum of squares, it suffices to restrict to the subspace 
{x = X\ = X2, y — Xs — 24, z = £5 = Xq} and show that the resulting polynomial 
W = (1 / 4) A 78 h(x , x, y, y, z, z, w, w) is not a sum of squares. This restriction is given by 

W = x 4 y 2 + 2x 3 y 3 + x 2 y 4 + x 4 yz + 5x 3 y 2 z + 6x 2 y 3 z + 2xy 4 z + x 4 z 2 + 5x 3 yz 2 + 10x 2 y 2 z 2 

+ 6xy 3 z 2 + y 4 z 2 + 2x 3 z 3 + 6x 2 yz 3 + Qxy 2 z 3 + 2y 3 z 3 + x 2 z 4 + 2xyz A + y 2 z 4 . 

This polynomial vanishes at six points in P 2 (IR), 

[1:0:0], [0:1:0], [0:0: 1], [1 : -1 : 0], [1:0: -1], and [0:1: -1]. 

Thus if W is written as a sum of squares h\, then each hk must vanish at each of these six 
points. The subspace of M.[x, y, z] 3 of cubics vanishing in these six points is four dimensional 
and spanned by v = {x 2 y + xy 2 , x 2 z + xz 2 , y 2 z + yz 2 , xyz}. Then W is a sum of squares if and 
only if there exists a positive semidefinite 4x4 matrix G such that W = v T Gv. However, 
the resulting linear equations in the variables Gy, 1 < % < j < 4, have the unique solution 

/ 1 1/2 1 2\ 
1/2 1 12 
1 112 
\ 2 2 2 5/ 

One can see that G is not positive semidefinite from its determinant, which is —1/4. Thus 
W cannot be written as a sum of squares. 

This, along with Corollary 4.3[ provides another proof that no power of the Vamos poly- 



G 



nomial h(x) has a definite determinantal representation. 

6. The cone of interlacers and its boundary 

Here we investigate the convex cone Int(/, e) of polynomials interlacing /. We compute 
this cone in two examples coming from optimization and discuss its algebraic boundary, the 
minimal polynomial vanishing on the boundary of Int(/, e), when this cone is full dimensional. 
For smooth polynomials, this algebraic boundary is irreducible. 

If the real variety of a hyperbolic polynomial / is smooth, then the cone Int(/, e) of 
interlacers is full dimensional in IRfx^-i- On the other hand, if V(/) has a real singular 
point, then every polynomial that interlaces / must pass through this point. This has two 
interesting consequences for the hyperbolic polynomials coming from linear programming 
and semidefinite programming. 

Example 6.1. Consider / = nr=i x «- ^ ne singular locus of V(f) consists of the set of 
vectors with two or more zero coordinates. The subspace of polynomials in IR[x] ra _i which 
vanish in these points is spanned by the n polynomials {flj^i x j '■ i = 1, • • • , fi}. Note that 



this is exactly the linear space spanned by the partial derivatives of /. Theorem |3.1| then 
shows that the cone of interlacers is isomorphic to C(f, e) = (R>o) n : 

int (n^i) = j x>n^ : ae ( R >°) n } - ( M >°)"' 

^ 1=1 j^i > 
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Interestingly, this also happens when we replace the positive orthant by the cone of positive 
definite symmetric matrices. 

Example 6.2. Let / = det(X) where X is a d x d symmetric matrix of variables. The 
singular locus of V(/) is the locus of matrices with rank < d— 2. The corresponding ideal is 
generated by the (d— 1) x (d— 1) minors of X. Since these have degree d— 1, we see that the 
polynomials interlacing det(X) must lie in the linear span of the (d — 1) x (d — 1) minors of 
X. Again, this is exactly the linear span of the directional derivatives De{J) = tr(E ■ X ad ->). 
Thus Theorem 3.1 identifies Int(/, e) with the cone of positive semidefinite matrices: 

Int (det(X), I) = I tr(A • X adj ) : A G Mg d j = 



!)dxd 
^0 • 



If Vr(/) is nonsingular, then the cone Int(/, e) is full dimensional and its algebraic bound- 
ary is a hypersurface in R[a;]d_i. We see that any polynomial g on the boundary of Int(/, e) 
must have a non-transverse intersection point with /. As we see in the next theorem, this 
algebraic condition exactly characterizes the algebraic boundary of Int(/, e). 

Theorem 6.3. Let f G R[x]d be hyperbolic with respect to e 6 R" and assume that the 
projective variety V(/) is smooth. Then the algebraic boundary of the convex cone Int(/, e) 
is the irreducible hypersurface in C[x]d-i given by 

(6.1) j# G C[x] d -i : 3 p G P n_1 such that f(p) = g(p) = and rank (y^jj < 1 



Proof. First, we show that the set (6.1) is irreducible. Consider the incidence variety X of 



polynomials g and points p satisfying this condition, 

X = E P(C[x] d _ 1 ) x P"- 1 : f(p) = g(p) = and rank (y^j) < 1 

The projection 7T2 onto the second factor is V(f) in P" --1 . Note that the fibres of ^ are linear 
spaces in P(C[x]d_i) of constant dimension. In particular, all fibres of iT2 are irreducible of 
the same dimension. Since X and V(/) are projective and the latter is irreducible, this 
implies that X is irreducible (see [17, §1.6, Thm. 8]), so its projection 7Ti(X) onto the first 



factor, which is our desired set (6.1), is also irreducible. 

If V(/) is smooth, then by |12| Lemma 2.4], / and D e f share no real roots. This shows 
that the set of polynomials g G K[x]d_i for which D e f ■ g is strictly positive on Vr(/) is 
nonempty, as it contains D e f itself. This set is open and contained in Int(/, e), so Int(/, e) 
is full dimensional in IR[x]d_i. Thus its algebraic boundary <9Int(/, e) is a hypersurface in 
C[x]rf_i. To finish the proof, we just need to show that this hypersurface is contained in 
( 6.1[ ), since the latter is irreducible 



To see this, suppose that g G lR[x]d-i lies in the boundary of Int(/, e). By Theorem 2.1 
there is some point p G Vr(/) at which g ■ D e f is zero. As / is nonsingular, D e f(p) cannot be 
zero, again using [T2| Lemma 2.4]. Thus g(p) = 0. Moreover, the polynomial g ■ D e f — f ■ D e g 
is globally nonnegative, so its gradient also vanishes at the point p. As f(p) = g(p) = 0, this 
means that D e f(p) ■ Vg(p) = D e g(p) ■ Vf(p). Thus the pair (g,p) belongs to X above. □ 
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When V(/) has real singularities, computing the dimension of Int(/, e) becomes more 
subtle. In particular, it depends on the type of singularity. 

Example 6.4. Consider the two hyperbolic quartic polynomials 

/i = 3y 4 + x 4 + 5x 3 z + 6x 2 z 2 — 6y 2 z 2 and f 2 = (x 2 + y 2 + 2xz)(x 2 + y 2 + 3xz), 

whose real varieties are shown in Figure [6] in the plane {z = 1}. Both are hyperbolic with 
respect to the point e = [— 1 : : 1] and singular at [0 : : 1]. Every polynomial interlacing 
either of these must pass through the point [0:0:1]. However, for a polynomial g to 
interlace / 2 , its partial derivative dg/dy must also vanish at [0 : : 1]. Thus Int(/i,e) has 
codimension one in M.[x,y, z\% whereas Int(/2, e) has codimension two. 




FIGURE 6. Two singular hyperbolic quartics with different dimensions of inter lacers. 



Theorem 



3.1 



states that C(f, e) is a linear slice of the cone Int(/, e). By taking boundaries 



of these cones, we recover V(/) as a linear slice of the algebraic boundary of Int(/, e). 

Definition 6.5. We say that a polynomial / G M.[x] is cylindrical if there exists a linear 
change of coordinates T on M. n such that f(Tx) G M[x\, . . . , x n _i]. 

Corollary 6.6. For non- cylindrical f , the map M. n — > M[x]d-i given by a t— >■ D a f is injective 
and maps the boundary ofC(f,e) into the boundary o/Int(/, e). If f is irreducible, this map 
identifies V(f) with a component of the Zariski closure of the boundary o/Int(/, e) in the 
plane spanned by df /dxi, . . . , df/dx n . 

Proof. Since / is not cylindrical, the n partial derivatives df/dxj are linearly independent, 
so that a i— > D a f is injective. The claim now follows from taking the boundaries of the cones 
in (3.1). If / is irreducible, then the Zariski closure of the boundary of C(f,e) is V(/). □ 

Example 6.7. We take the cubic form f(x, y, z) = {x — y)(x + y)(x + 2y) —xz 2 , which is hy- 
perbolic with respect to the point [1:0:0]. Using the computer algebra system Macaulay 2, 
we can calculate the minimal polynomial in Q[cn, C\2, C13, C22, C23, C33] that vanishes on the 
boundary of the cone Int(/, e). Note that any conic 

q = cnx 2 + c 12 xy + c 13 xz + c 2 iy 2 + c 23 yz + c 33 z 2 

on the boundary of Int(/, e) must have a singular intersection point with V(/). Saturating 
with the ideal (x, y, z) and eliminating the variables x, y, z from the ideal 

(/) <?) + minors2(Jacobian(/, q)) 
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gives an irreducible degree-twelve polynomial in the six coefficients of q. This hypersurface 
is the algebraic boundary of Int(/, e). When we restrict to the three-dimensional subspace 
given by q — + + c|£, this degree-twelve polynomial factors as 

a ■ f(a, 6, c) • (961a 8 + 5952a 7 6 + 11076a 6 6 2 - 3416a 5 6 3 - 34770a 4 6 4 - 31344a 3 6 5 
+ 14884a 2 6 6 + 34632a6 7 + 136 8 96 s - 1896a 6 c 2 - 4440a 5 6c 2 + 6984aW + 25728a 3 6 3 c 2 
+ 15960aW - 7560a6 5 c 2 - 75606 6 c 2 + 1074a 4 c 4 - 1680a 3 6c 4 - 7116a 2 6 2 c 4 - 2376a6 3 c 4 
+ 21066V + 16a 2 c 6 + 936a6c 6 - 27c 8 ). 

One might hope that Int(/, e) is also a hyperbolicity cone of some hyperbolic polynomial, 
but we see that this is not the case. Restricting to c = shows that the polynomial above, 
unlike /, is not hyperbolic with respect to [1 : : 0]. 
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